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ABSTRACT 

For laminated ceramic matrix composite (CMC) mate- 
rials to realize their full potential in aerospace 
applications, design methods and protocols are a nec- 
essity. This paper focuses on the t ime- independent 
failure response of these materials and presents a 
reliability analysis associated with the initiation of 
matrix cracking. It highlights a public domain computer 
algorithm that has been coupled with the laminate analy- 
sis of a finite element code and which serves as a 
design aid to analyze structural components made from 
laminated CMC materials. Issues relevant to the effect 
of the size of the component are discussed, and a param- 
eter estimation procedure is presented. The estimation 
procedure allows three parameters to be calculated from 
a failure population that has an underlying Weibull 
distribution. 

INTRODUCTION 

Structural components produced from laminated 
ceramic matrix composite (CMC) materials are being con- 
sidered for a broad range of aerospace applications that 
include propulsion subsystems in the national aerospace 
plane, the space shuttle main engine, and advanced gas 
turbine engines. Specifically, composite ceramics may 
be used as segmented engine liners, small missile engine 
turbine rotors, and exhaust nozzles. These materials 
will improve fuel efficiency by increasing engine tem- 
peratures and pressures, which will, in turn, generate 
more power and thrust. Furthermore, these materials 
have significant potential for raising the thrust-to- 
weight ratio of gas turbine engines by tailoring 
directions of high specific reliability. The emerging 
composite systems, particularly those with a silicon 
nitride or silicon carbide matrix, can compete with 
metals in many demanding applications. The capabili- 
ties of laminated CMC prototypes have already been 
demonstrated at temperatures approaching 1400 °C, well 
beyond the operational limits of most metallic 
materials . 
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Adding a second ceramic phase with an optimized 
interface to a brittle matrix improves fracture tough- 
ness, decreases the sensitivity of the brittle ceramic 
matrix to microscopic flaws, and could also improve 
strength. The presence of fibers in the vicinity of 
the crack tip modifies fracture behavior by increasing 
the required crack driving force by several mechanisms. 
These mechanisms include crack pinning, fiber bridging, 
fiber debonding, and fiber pull-out. This increase in 
fracture toughness allows for ^graceful* rather than 
catastrophic failure. A unidirectional ply loaded in 
the fiber direction retains substantial strength capac- 
ity beyond the initiation of matrix cracking despite the 
fact that neither of the constituents would exhibit such 
behavior if tested alone. First matrix cracking consis- 
tently occurs at strains greater than in the monolithic 
matrix material. As additional load is applied, the 
matrix tends to break in a series of cracks bridged by 
the ceramic fibers, until the ultimate strength of the 
composite is reached. The region of a typical stress- 
strain curve between the first matrix cracking and the 
ultimate tensi le strength illustrates an intrinsic dam- 
age tolerance not present in monolithic ceramics. 

Laminated CMC material systems have several mechan- 
ical characteristics which must be considered in the 
design of structural components. In regards to an indi- 
vidual ply, the most deleterious of these characteris- 
tics are low strain tolerance, low fracture toughness, 
and a large variation in failure strength in the mate- 
rial orientation transverse to the fiber direction. 

Thus analyses of components fabricated from ceramic 
materials require a departure from the usual determinis- 
tic philosophy of designing metallic structural compo- 
nents (i.e., the f actor-of-saf ety approach). Although 
the so-called size effect has been reported to be non- 
existent in the fiber direction (see Dicarlo, 1989), the 
bulk strength of unidi rect ional-reinforced ply will 
decrease transverse to the fiber direction as the compo- 
nent volume increases. Since failure in the transverse 
direction will be dominated by the scatter in strength, 
statistical design approaches must be employed. These 
approaches must, on the one hand, allow for elevated 
strength, reduced variability in strength, and a dimin- 
ished effect of bulk specimen size in the fiber direc- 
tion, and, on the other hand, increased scatter in 
strength and effects of bulk size in the transverse 


direction. Simply stated, a reliability analysis must 
rationally account for material symmetry imposed by the 
reinforcement. Computational structural mechanics phil- 
osophies must emerge that address the issues of scatter 
in strength, size effect, and material anisotropy. 

There is a need for test-bed software programs that 
incorporate stochastic design protocols, that are user 
friendly, that are computationally efficient, and that 
have flexible architectures which can readily incorpo- 
rate changes in design philosophy. The C/CARES (Compos- 
ite Ceramics Analysis and Reliability Evaluation of 
Structures) program, which will be highlighted in this 
article, was developed to fulfill this need. C/CARES is 
a public domain computer algorithm, coupled to a general 
purpose finite element program, which predicts the fast 
fracture reliability of a laminated structural component 
under multiaxial loading conditions. 

STOCHASTIC DESIGN ISSUES 

For a number of composite material systems, sev- 
eral authors (see for example Batdorf and Chaffarian, 
1984; and Wu, 1989) have reported a diminished size 
effect in the fiber direction; and DiCarlo (1989) has 
reported this effect for ceramic composites, in parti- 
cular. This phenomenon is an important feature that 
must be addressed by any reliability model. How it is 
addressed depends on whether the material is modeled as 
a series system, a parallel system, or a combination. 
Current analytical practice uses finite element methods 
to determine the state of stress throughout the compo- 
nent. It is assumed that failure depends on the stress 
state in a component, such that deformations are not 
controlling design. Since failure may initiate in any 
of the discrete volumes (elements), it is useful to con- 
sider a component from a systems viewpoint. A component 
comprised of discrete volumes is a series system if it 
fails when one of the discrete volumes fail. This 
approach gives rise to weakest-link theories. In a 
parallel system, failure of a single element does not 
necessarily cause the component to fail since the 
remaining elements may be able to sustain the load 
through redistribution. The parallel system approach 
leads to what has been referred to in the literature as 
"bundle” theories. 

The basic principles underlying these bundle theo- 
ries were originally discussed by Daniels (1945) and 
Coleman (1958). Their work was extended to polymer 
matrix composites by Rosen (1964) and Zweben (1968). 
Here, a relatively soft matrix serves to transfer stress 
between fibers and contributes little to the composite 
tensile strength. Hence, when a fiber breaks the load 
is transferred only to neighboring fibers. Their analy- 
sis is rather complex and limited to establishing bounds 
on the stress at which the first fiber breaks and the 
stress at which all the fibers are broken. Harlow and 
Phoenix (1978) proposed a rather abstract approach that 
established a closed-form solution for all the interme- 
diate stress levels in a two-dimensional problem, and 
Batdorf (1982) used an approximate solution to estab- 
lish the solutions for the three-dimensional problem. 
Batdorf’ s model includes the two-dimensional model as a 
special case. In both of the latter two models, the 
authors proposed that the effective Weibull modulus 
increases with increasing component volume. This 
implies a diminished size effect. However, these cur- 
rent bundle theories are predicated on the fact that 
fibers are inherently much stronger and stiffer than the 
matrix. In laminated CMC materials this is not always 
the case. The strength and stiffness of both the fiber 
and matrix are usually closer in magnitude. For this 
reason bundle theories will not be considered in this 
paper . 


We advocate the use of a weakest-link reliability 
theory for designing components manufactured from lami- 
nated CMC materials that do not exhibit strong size 
effects in specific directions. Assuming that a lami- 
nated structure behaves in a weakest-link manner allows 
a conservative estimate of structural reliability to be 
calculated. Thomas and Wetherhold (1990) point out that 
this assumption is equivalent to predicting the proba- 
bility of the first matrix crack occurring in an indi- 
vidual ply. For most applications the design failure 
stress for a laminated material is assumed to coincide 
with this first ply matrix cracking because matrix 
cracking usually allows the fibers to oxidize at high 
temperatures, embrittling the composite. 

Next, we address the righteousness of applying 
weakest-link theory to a material that in some sense 
does not exhibit size effects. In general, the mean 
strength of a sample population representing uniaxial 
tension test specimens can be obtained by integrating 
the probability of survival P s with respect to the 
applied tensile stress: that is, 

a = [ P da ( 1 ) 

J 0 

Here a is the mean tensile strength, a is the applied 
tensile stress, and P s is the probability of survival, 
The form of P s depends on the probability density 
function that best represents the failure data, which in 
turn depends on whether the structural component acts as 
a parallel or series system. Adopting a three-parameter 
Weibull cumulative distribution function, and assuming 
a weakest-link system (a conservative assumption) gives 
the following form to Eq. (1): 
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Here V is the volume of the tensile test specimen. 

(3 is the scale parameter, a is the Weibull shape 
parameter, and y is the threshold stress. This inte- - 
gral has the following closed-form solution (see 
DeSalvo 1970): 



which depends on volume, the Weibull parameters, and the 
gamma function T. When an argument originally outlined 
by Jayatilaka (1979) is followed, two uniaxial tensile 
specimen populations with distinctly different specimen 
volumes will yield different mean strengths. Associat- 
ing a i with gives 



Similarly associating 02 w ith V 2 gives 



If the effective mean is defined as 
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then the ratio of the effective mean strengths depends 
only on the specimen volume and the Weibull modulus: 
that is, 



As the Weibull modulus of a particular material 
increases, the ratio of the effective mean strengths 
approaches unity. In this situation the material 
exhibits no size effect (even though the distribution 
of failure strength may be represented by a Weibull 
probability density function). From a practical stand- 
point, doubling the specimen size of a material whose 
Weibull modulus is =15 would yield less than a 5 per- 
cent difference in the effective mean failure strengths 
of the two populations. We expect an elevated Weibull 
modulus to be associated with the strength of CMC mate- 
rials in the fiber direction. Reports of an apparent 
lack of size effect associated with the strength in the 
fiber direction (see DiCarlo, 1989) could easily be an 
artifact of an increasing shape parameter (or small 
sample size). However, at this time there is an 
insufficient quantity of CMC failure data from which 
to estimate the Weibull parameters. In general, the 
weakest-link theory allows for diminishing size effects 
as the Weibull modulus increases. 

RELIABILITY MODEL 

The ongoing metamorphosis of ceramic material sys- 
tems and the lack of standardized design data has in the 
past tended to minimize the emphasis on modeling. Many 
structural components fabricated from ceramic materials 
were designed by "trial and error,” since emphasis was 
placed on demonstrating feasibility rather than on fully 
understanding the processes controlling behavior. (This 
is understandable during periods of rapid improvements 
in material properties for any system.) In predicting 
failure behavior, there is a philosophical division that 
separates analytical schools of thought into microstruc- 
tural methods (usually based on principles of fracture 
mechanics) and phenomenological methods. Blass and 
Ruggles (1990) point out that analysts from the first 
school would design the material assuming that the con- 
stituents are distinct structural components and would 
consider the composite ply (or lamina) a structure in 
its own right. Analysts from the latter school of 
thought would design with the material (i.e., they would 
analyze structural components fabricated from the mate- 
rial). Rigorous fracture mechanics criteria have been 
proposed (e.g., Budiansky et al., 1986; and Marshall 
et al . , 1985) that adopt the microstructural viewpoint, 
but since they are all deterministic criteria, they will 
not be considered here. Fracture mechanics has been 
combined with a probabilistic Weibull analysis of fail- 
ure location to determine the stress-strain behavior and 
subsequent work of fracture for unidirectional compos- 
ites (e.g., Thouless and Evans, 1988; and Sutcu, 1989). 
However, the focus here is first matrix cracking, and 
we note that mature reliability-based design methods 
using fracture mechanics concepts will not surface un- 
til a coherent mixed-mode fracture criterion has been 
proposed. 

The aforementioned second school of thought repre- 
sents the ply (or lamina) as a homogenized material 
with strength properties that are determined from a 


number of well-planned phenomenological experiments. 

The authors currently embrace this philosophy, and there 
are practical reasons for initially adopting this view- 
point. We fully recognize that the failure characteris- 
tics of these composites are controlled by a number of 
local phenomena including matrix cracking, debonding 
and slipping between matrix and fibers, and fiber break- 
age, all of which strongly interact. Understanding the 
underlying analytical concepts associated with the 
microstructural viewpoint allows one to gain insight 
and intuition prior to constructing muitiaxial failure 
theories that in some respect reflect the local behav- 
ior. Tensile failure in the fiber direction is depend- 
ent on these local mechanisms, and the future intent is 
to extend reliability methods to the constituent level 
in a rational and practical manner. However, a top-down 
approach, that is first proposing design models at the 
ply level, will establish viable and working design pro- 
tocols. Initially adopting the bottom-up approach 
allows for the possibility of becoming mired in detail 
(experimental and analytical) when muitiaxial reliabil- 
ity analyses are conducted at the constituent level. 

There is a great deal of intrinsic variability in 
the strength of each brittle constituent of a ceramic 
matrix composite, but depending on the composite system, 
the transverse matrix cracking strength may either be 
deterministic or probabilistic. Statistical models are 
a necessity for those composite systems which exhibit 
any scatter in the initiation of first matrix cracking. 
We treat it in a probabilistic fashion, requiring that 
deterministic strength be a limiting case that is read- 
ily obtainable from the proposed reliability model. 
Predicting the reduction in reliability due to loads in 
the fiber direction addresses an upper bound for ply 
reliability in a structural design problem. Conversely, 
a tensile load applied transverse to the fiber direction 
results in failure behavior similar to a monolithic 
ceramic, which corresponds to the lower bound of ply 
reliability. Thus muitiaxial design methods must be 
capable of predict ing Jhese two bounds as well as 
account for the reduction in reliability due to an 
in-plane shear stress, and compressive stresses in the 
fiber direction and transverse to the fiber direction. 

A number of macroscopic theories exist that treat 
unidirectional composites as homogenized, anisotropic 
materials. These methods use phenomenological strength 
data directly without hypothesizing specific crack 
shapes or distributions. Theories of this genre gener- 
ally are termed noninteractive if individual stress com- 
ponents are compared to their strengths separately. In 
essence, failure mechanisms are assumed not to interact, 
and this results in component reliability computations 
that are quite tractable. Work by Thomas and Wetherhold 
(1990), Duffy and Arnold (1990), Duffy and Manderscheid 
(1990), and Duffy et al . (1990) are representative of 
muitiaxial noninteractive reliability models for aniso- 
tropic materials. In addition Wu (1989), and Hu and 
Coetschel (1989) have proposed simpler unidirectional 
reliability models for laminated composites that can be 
classified as noninteractive. Alternatively, one can 
assume that for muitiaxial states of stress, failure 
mechanisms interact and depend on specific stochastic 
combinations of material strengths. Usually a failure 
criterion is adapted from exis t i ng polymer matrix design 
technologies. The probability that the criterion has 
been violated for a given stress state is computed using 
Monte Carlo methods (de Roo and Paluch, 1985) or first- 
order-second-moment (F0SM) methods (Yang, 1989, and Miki 
et al., 1989). The interactive approach often results 
in computationally intensive reliability predictions. 

In this paper a noninteractive phenomenological 
approach has been chosen such that a unidirectional ply 
is considered a two-dimensional structure, assumed to 
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have five basic strengths (or failure modes). They 
include a tensile and compressive strength in the fiber 
direction, a tensile and compressive strength in the 
direction transverse to the fiber direction, and an 
in-plane shear strength. In addition each ply is dis- 
cretized into individual sub-ply volumes. For reasons 
discussed in the previous section we assume that fail- 
ure of a ply is governed by its weakest link (or sub-ply 
volume). Under this assumption, events leading to fail- 
ure of a given link do not affect other links (see, for 
example, Batdorf and Heinisch, 1978; Wetherhold, 1983; 
and Cassenti, 1984); thus the reliability of the ith ply 
is given by the following expression 



where V is the component volume. Here, ^i(xj) is the 
failure function per unit volume at position xj within 
the ply , given by 
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The a's associated with each term in Eq. (9) correspond 
to the Weibull shape parameters, the fi's correspond to 
Weibull scale parameters, and the y s correspond to the 
Weibull threshold stresses. In addition, a* and <73 
represent the in-plane normal stresses that are aligned 
with and transverse to the fiber direction, respec- 
tively. Also, x i 2 IS in-plane shear stress. The 
normal stresses appear twice and this allows for differ- 
ent failure modes to emerge in tension and compression. 
Mote that the brackets indicate a unit step function: 
i .e. , 

( x x > 0 

(x> = x • u[x] = < ( 10 ) 

(0 x < 0 


87.5 MPa were applied to the tube. The Weibull param- 
eters were also arbitrarily chosen (see Table 1 ). Mote 
that the threshold stresses are taken as zero for sim- 
plicity. In design, setting the threshold stresses 
equal to zero would represent a conservative assumption. 
The overall component reliability is depicted as a func- 
tion of the midply orientation angle (0) in Fig. 2. 

The ply orientation has a decided effect on component 
reliability, as expected. Similar studies could demon- 
strate the effects of component geometry, ply thick- 
ness, load, and/or Weibull parameters on component reli- 
ability. Hence, the C/CARES code allows the design 
engineer a wide latitude to optimize a component rela- 
tive to a number of design parameters. 

PARAMETER ESTIMATION 

We anticipate that laminated CMC materials will 
exhibit threshold behavior, at least in the fiber direc- 
tion. Hence, a three-parameter Weibull distribution is 
used in the stochastic failure analysis of the compo- 
nents. The threshold stress parameter is included to 
allow for zero probability of failure when the load is 
below a predetermined level. The three-parameter dis- 
tribution has been somewhat ignored due to difficulties 
encountered in extracting the parameters from experi- 
mental data. Several authors (including Weibull, 1939; 
Weil and Daniel. 1964; and Schneider and Palazotto. 

1979) have proposed estimation methods for the three- 
parameter distribution. For various reasons, these 
techniques have not been widely accepted. However, 
Cooper (1988) recently proposed a nonlinear regression 
method to estimate parameters. Regression analysis pos- 
tulates a relationship between two variables. In an 
experiment, typically one variable can be controlled 
(the independent variable) while the response variable 
(or dependent variable) is uncontrolled. In simple 
failure experiments the material dictates the strength 
at failure, indicating that the failure stress is the 
response variable. The ranked probability of failure 
* i can be controlled by the experimentalist since it 
is functionally dependent on the sample size M. If 
the observed failure stresses (a^, < 73 , < 73 , * . . , <t^) 
are placed in ascending order, then 

W - irfti (12) 


Inserting Eq. ( 9 ) into the volume integration given by 
Eq. ( 8 ) yields the reliability of the ith ply, and the 
probability of first ply failure for the laminate is 
given by the expression 


p fpf 
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where n is the number of plies. 

This reliability model can be readily integrated 
with laminate analysis options available in several com- 
mercial finite element codes. A preliminary version of 
a public domain computer algorithm (C/CARES) that is 
coupled with MSC/NASTRAN has been developed at NASA 
Lewis Research Center to perform this analysis. A sim- 
ple benchmark application illustrates the approach. A 
thin-wall tube is subjected to an internal pressure and 
an axial compressive load. The component is fabricated 
from a three-ply laminate, with a 9O°/0/9O° layup. Here 
angle 0 is measured relative to the longitudinal axis 
of the tube (see Fig. 1). An arbitrary internal pres- 
sure of 4.25 MPa and an axial compressive stress of 


Clearly one can influence the ranked probability for a 
given stress level by increasing or decreasing the sam- 
ple size. Cooper's procedure adopts this philosophy, 
and the specimen failure stress is treated as the 
dependent variable. The associated ranked probability 
of failure is assumed to be the independent variable. 
The basic three-parameter Weibull expression for proba- 
bility of failure can be expressed as 
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where a. is an estimate of the dependent variable, and 

Y, D, and a are estimates of the threshold parameter, 
the characteristic strength, and the shape parameter, 
respectively. Defining 
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as the ith residual, where as before a. is the ith 
failure stress, then 
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where we adopt Cooper's notation and take 

Setting the partial derivatives of the sum of the 

squares of the residuals with respect to y, 3, and a 
equal to zero yields the following three expressions: 
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in terms of the parameter estimates. The solution of 
this system of equations is iterative. One assumes an 
initial value for a (a small value, usually equal 

to 1), computes 3 from Eq, (17) and y from Eq. (18). 
These values of the parameter estimates are then inser- 
ted into Eq. (19), and this expression is checked to 
see if it satisfies some predetermined tolerance. If 
Eq. (19) is not satisfied, a is increased and a new 
iteration is conducted. This procedure continues until 
a set of parameter estimates are determined that satisfy 
Eqs . (17) to (19). 

Currently we do not have enough CMC failure data to 
estimate Weibull parameters for a given material orien- 
tation. So to illustrate the technique, parameter 
estimates were determined for two-parameter and three- 
parameter distributions from a failure population repre- 
senting a monolithic ceramic (alumina) reported by Quinn 
(1989). The failure data and specimen geometry are 
shown in Table 2 and Fig. 3. Figure 4 is a plot of 
probability of failure versus failure stress for the 
data. The straight line represents the two-parameter 
fit to the data, using Quinn’s (1989) values for the 
shape and scale parameters. The nonlinear curve repre- 
sents the three-parameter fit to the data. Note that 
the three-parameter distribution is more efficient in 


predicting the failure data in the high reliability 
region. 

SUMMARY AND FUTURE DIRECTIONS 

In this paper we discuss stochastic issues related 
to size effects in the fiber direction of a unidirec- 
tional CMC material. In addition, we present a relia- 
bility model along with a simple application that high- 
lighted the C/CARES computer algorithm. (This public 
domain algorithm is capable of predicting component 
reliability from the state of stress and temperature 
distribution within the component.) The authors antic- 
ipate that CMC materials will exhibit threshold behav- 
ior, hence a nonlinear regression analysis was outlined 
to determine three parameters for a Weibull distribu- 
tion from failure data. 

Ceramic material systems will play a significant 
role in future elevated-temperature applications. To 
this end, there are a number of issues that must be 
addressed by the structural mechanics research commu- 
nity. We begin by pointing out that total failure of an 
individual ply effectively reduces the overall laminate 
stiffness. This causes local redistribution of the load 
to adjacent layers. In addition, delamination between 
plies relaxes the constraining effects among layers, 
allowing in-plane strains to vary in a stepwise fashion 
within a laminate. These effects require the develop- 
ment of rational load redistribution schemes. It is 
also apparent that before ceramics are used as struc- 
tural components in harsh service environments, thought- 
ful consideration should be given to reliability 
degradation due to time-dependent phenomena. Thus, 
issues germane to component life, such as cyclic fatigue 
and creep behavior, must be addressed analytically. 
Computational strategies are needed to extend current 
methods of analysis from subcritical crack growth and 
creep rupture to laminated CMC materials that are sub- 
ject to multiaxial states of stress. 

An important aspec_t that has not been addressed in 
detail is the effect of a rising R-curve behavior, where 
fracture toughness is a function of crack size. Clearly 
f iber- toughened matrices have process zones around the 
crack tip. Within this zone, energy dissipates locally, 
which develops a damage tolerance by increasing the 
resistance to crack growth with crack extension. Fail- 
ure of materials exhibiting R-curve behavior would 
depend on the rate at which resistance increases with 
crack growth. During crack extension this behavior 
would modify the strength distribution. Modeling 
efforts by Kendall et al. (1986) and others have 
accounted for this behavior in monolithic ceramics, and 
it is reasonable to extend their work to ceramic compos- 
ite material systems. Furthermore, if ceramic materials 
mimic ductile failure locally, cyclic fatigue may become 
a design issue, tinder cyclic loads, the process zone 
advances as the crack tip extends; therefore, brittle 
fracture mechanics may need to be modified to account 
for pseudoduc t i le fracture. Hence, application of 
modified metallic fatigue analyses may be a distinct 
possibility. 

In addition, recent progress in processing ceramic 
material systems has not been matched by mechanical 
testing efforts. There is a definite need for experi- 
ments that support the development of reliability mod- 
els. Initially this effort should include experiments 
that test the fundamental concepts (e.g., quantifying 
the size effect in the fiber direction) within the 
framework of current stochastic models. For example, 
probing experiments could be conducted along various 
biaxial load paths to establish level surfaces of reli- 
ability in a particular two-dimensional stress space 
(similar to probing yield surfaces in metals). One 
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could then verify such concepts as the maximum stress 
response which is often assumed in the noninteractive 
reliability models proposed for these materials. After 
establishing a theoretical framework, characterization 
tests should then be conducted to provide the functional 
dependence of model parameters with respect to tempera- 
ture and environment. Finally data from structural 
tests that are multiaxial (and possibly noni so thermal ) 
would be used to challenge the predictive capabilities 
of models through comparison to benchmark response data. 
These tests involve nonhomogeneous fields of stress, 
deformation, and temperature, and would include two-bar 
tests as well as plate and shell structures. Results 
from structural testing provide feedback for subsequent 
modification, but ad hoc models result in the absence of 
structured interaction between the experimentalist and 
the theoretician. The validity of these models are then 
forever open to question. Furthermore, we cannot over- 
emphasize that this kind of testing supports methods for 
designing components, not the materials. Currently this 
effort is hampered by the quality and scarcity of data 
(note the lack of failure data necessary to estimate 
composite Weibull parameters). Finally, ceramic pro- 
perties pertinent to structural design (which include 
stochastic parameters) vary with test methods. The 
mechanics research community is beginning to realize 
this, and a consensus is beginning to form regarding 
standards. However, we wish to underscore the funda- 
mental need for experimental programs that are relevant 
to structural mechanics issues. 

In closing, we recognize that when failure is less 
sensitive to imperfections in the material, stochastic 
methods may not be as essential. Yet, trends in design 
protocols are moving in the direction of probabilistic 
analyses (even for metals) and away from the simplistic 
safety-factor approach. In this sense, brittle ceramics 
will serve as prototypical materials in the study and 
development of reliability models that will act as the 
basis of future design codes. 
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TABLE 1. - COMPOSITE WEIBULL PARAMETERS FOR THIN-WALL TUBE 


[Weibull threshold stress, y[ - 0.] 


Index a 

Type and direction of stress 

Weibull parameters 



Shape . 
a i 

Scale , 
Gi 

1 

Normal tensile stress in fiber direction 

25 

450 

2 

In-plane shear stress 

22 

420 

3 

Normal tensile stress transverse to 
fiber direction 

to 

350 

4 

Normal compressive stress in fiber 
direction 

35 

4500 

5 

Normal compressive stress transverse to 
fiber direction 

30 

3500 


a Indices correspond to subscripts in Eq. (9). 



Fig. 1 . Finite element mesh of thin-wal! tube with three-ply 
{90°/6/90°) layup. Internal pressure, 4.25 MPa; axial com- 
pressive stress, 87.5 MPa. 
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Overall component reliability 


1.0 



Midply orientation angle 

Fig. 2. Component reliability versus midply orientation 
angle for thin wall tube. 



20 mm 



Number 

of 

Shape, 

Weibuil parameters 
Scale, Threshold stress, 

parameters 

a 

P 

V 

3 

1.15 

803.41 

298.48 

2 

13.2 

376.0 



o Data (35 points) 



Failure stress, MPa 

Fig. 4. Two-parameter and three-parameter distrib- 
utions determined from the alumina failure. 


TABLE 2. - MONOLITHIC ALUMINA FAILURE DATA 3 


Speci- 

men 

number 

Stress , 
MPa 

28 

385 

29 

388 

30 

395 

31 

402 

32 

411 

33 

413 

34 

415 

35 

456 


Speci- 

men 

number 

Stress , 
MPa 


Speci- 

men 

number 

Stress , 
MPa 

10 

337 


19 

357 

11 

343 


20 

364 

12 

3 45 


21 

371 

13 

347 


22 

373 

14 

350 


23 

374 

15 

352 


24 

375 

16 

353 


25 

376 

17 

355 


26 

376 

18 

356 


27 

381 


Speci- 

men 

number 

Stress , 
MPa 

1 

307 

2 

308 

3 

322 

4 

328 

5 

328 

6 

329 

7 

331 

8 

332 

9 

335 


a For specimen shown in Fig. 3. 
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